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Abstract 

Using data obtained with the CLEO III detector, running at the Cornell Electron Storage Ring 
(CESR), we report on a new study of exclusive radiative T(1S) decays into the final states 77r + 7r~, 
jK + K~, and 7pp. We present branching ratio measurements for the decay modes T(1S) — > 
7/2(1270), T(1S) -► 7/^(1525), and T(1S) -► jK+K'; helicity production ratios for / 2 (1270) and 
/^(1525); upper limits for the decay T(1S) -► 7/j(2200), with /j(2220) -► ir + ir~, K+K~, pp; and 
an upper limit for the decay Y(1S) — > 7^(1860), with X(1860) — > 7pp. 
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I. INTRODUCTION 



Radiative decays of heavy quarkonia, where a photon replaces one of the three gluons 
from the strong decay of, for example, the J/ip or T(1S), are useful in studying color-singlet 
two-gluon systems. The two gluons can, among other things, hadronize into a meson 1 , or 
directly form a glueball 2 . Further information on radiative decays of heavy quarkonia can 
be found in Q . 

Light-meson production in J /if) two-body radiative decays has been experimentally well 
established with branching fractions at the 10~ 3 level, based largely on evidence provided by 
radiative decays to a pair of hadrons 3 . The production ratios of the available helicity states 



have been measured for the tensor mesons ^(1270) [46J, |47|, |48j, |49j and f2^1525) [5CJ, |5l| 
in J /if) two-body radiative decays and agree with theoretical predictions [52[[53j]. In 1996, 
the BES Collaboration reported the observation of the /j(2220) in J /if) two-body radiative 
decays, and measured product branching fractions, B(J/ip — > 7/j(2220)) x £>(/j(2220) — > 
h + h~) (we use the convention h = ir, K, p), of the order of 10 -5 jHiJ]. Much excitement 
was generated at the time because it is possible to interpret the /j(2220) as a glueball. A 
candidate similar to /j(2220) was reported in 1986 by the Mark III collaboration in the KK 
mode [HH, but was not confirmed by the DM2 collaboration [56|. Recently, BES reported the 
existence of a new particle, the X(1860), observed by its d ecay J / if) — » 7X(1860) — > pp |5?| . 
a result that is currently being interpreted [H,[5^, EkJ 61, 62[ 63L 164^. 

The experimental observation of radiative T(1S) decays is challenging because their rate 
is suppressed to a level of 

'qt V ( rric 
Ac) \m K 

of the corresponding rate of J/if> radiative decays. This factor arises because the quark- 
photon coupling is proportional to the electric charge, and the quark propagator is rou ghly 
proportional to 1/m for low momentum quarks. Taking into account the total widths [45j 
of J /if) and T(1S), the branching fraction of a particular T(1S) radiative decay mode is 
expected to be around 0.04 of the corresponding Jif) branching fraction. In 1999, CLEO II 
made the first observation of a radiative T(1S) decay to a pair of hadrons 65], which 
was consistent with T(1S) — > 7/2(1270), where ^(1270) — > tttt. Comparing the measured 
branching fraction to the J /if) — > 7/2(1270) branching fraction, a suppression factor of 



0.025, 



1 Several authors have studied meson production in Y(1S) radiative decays, giving predictions for branching 
and helicity production ratios. The heavy-quarkonium system is usually described by non-relativistic 
QCD while the gluonic hadronization has been treated using soft collinear effective theory gluon 
distribution amplitudes and perturbative QCD 0,0). 

2 Glueballs are a natural consequence of QCD, and predictions of their properties have been made using 
different approaches, such as potential models BH1, lattice QCD calculations HHQEl 

, bag mod- 
els El El 0], flux-tube models E3, the QCD sum rules Q, the Bethe-Salpeter (B-S) equation 
, QCD factorization formalism models [^llE^I . weakly-bound-state models 01, and a three-dimensional 
relativistic equation 0. However, despite intense experimental searches [25L |2qL E^i Esi El 0. lilf . there 
is no conclusive experimental evidence of their direct observation, alt houg h there are strong indications 
that glueballs contribute to the rich light scalar and tensor 42, 431 

spectrums. 

3 We refer to the Particle Data Group [45j for a summary of J/ip radiative decays. 
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0.06 ±0.03 was obtained. Recent theoretical work [2j,|3j] predict a suppression factor between 
0.06 — 0.18 for this mode, and favor the production of /2(1270) in a helicity-0 state. After the 
BES result for the /j(2220) in radiative J/ib decays, a corresponding search was performed 
by CLEO II in the radiative T(1S) system [66j and limits were put on some of the glueball 
candidates' product branching ratios. 

In this paper, we use the CLEO III T(1S) data sample, which has fifteen times higher 
statistics and better particle identification than the CLEO II data sample, to probe the 
color-singlet two-gluon spectrum by measuring the system's invariant mass using its decays 
to 7r + 7r~, K + K~ , and pp. Further details of this analysis can be found elsewhere |67| . 



II. CLEO III DETECTOR, DATA, AND MONTE CARLO SIMULATED SAMPLE 

The CLEO III detector is a versatile multi-purpose particle detector described more fully 
in (68[. It is centered on the interaction region of CESR. From the e + e~ interaction region 
radially outward it consists of a silicon strip vertex detector and a wire drift chamber used 
to measure the position, momenta, and specific ionization energy losses (dE/dx) of charged 
tracks based on their fitted path in a 1.5 T solenoidal magnetic field and the amount of charge 
deposited on the drift chamber wires. The silicon vertex detector and drift chamber tracking 
system achieves a charged particle momentum resolution of 0.35%(1%) at 1 GeV/c(5 GeV/c) 
and a dE/dx resolution of 6%. Beyond the drift chamber is a Ring Imaging Cherenkov 
Detector, RICH, which covers 80% of the solid angle and is used to further identify charged 
particles by giving for each mass hypothesis the likelihood of a fit to the Cherenkov radiation 
pattern. After the RICH is a Crystal Calorimeter (CC) that covers 93% of the solid angle. 
The CC has a resolution of 2.2% (1.5%) for lGeV(5GeV) photons. After the CC is a 
superconducting solenoid coil that provides the magnetic field, followed by iron flux return 
plates with wire chambers interspersed in three layers at 3, 5, and 7 hadronic interaction 
lengths to provide muon identification. 

The data sample has an integrated luminosity of 1.13 fb 1 taken at the T(1S) energy, 
\fs = 9.46 GeV, which correspond to 21.2±0.2 million T(1S) decays [69] and 3.49 fb -1 taken 
at the T(4S) energy, y/s = 10.56 GeV, used to model the underlying continuum present in the 
T(1S) data sample. The continuum background modeling is important because continuum 
background processes such as e + e~ —>■ 7p with p — > 7r + 7r~, e + e~ — > 70 with — > K + K~ , 
and direct e + e~ — ► 'yh + h~ have the same topology as the signal events we are investigating. 

Efficiencies are evaluated using a Monte Carlo simulation of the process [7(J and a 
GEANT-based jlH detector response. Monte Carlo samples of e + e~ — > 7X with X — > h + h~ 
are generated at both the T(1S) and T(4S) energies with uniform angular distributions and 
flat h + h~ invariant mass distributions from threshold to 3.5GeV/c 2 . 



III. EVENT SELECTION 

Events which satisfy the CLEO III trigger [z^ are then required to meet the following 
analysis requirements: (a) There are exactly two charged tracks that trace back to the 
beamspot and have good quality track fits and dE/dx information, (b) There is exactly one 
CC shower that is unmatched to any track and whose energy, E~, is greater than 4 GeV. 

Each event is also required to be consistent with having the 4-momentum of the ini- 
tial e + e~ system by demanding that the chi-squared from a kinematic fit to the following 
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constraint, 

Ph+h- + (2E heiim — E h + h -)p^ = pom, (1) 

be less than 100, where Ph+h- is the di-hadron momentum, E h + h - is the di-hadron energy, 
Ebeum is the beam energy, p 7 is the photon's direction, and pcm is the momentum of the e + e~ 
system (which has a magnitude of a few MeV/c because of the small (~ 2 mrad) crossing 
angle of the e + and e~ beams). Equation Q is a 3-constraint subset of the 4-momentum 
constraint and has the convenient property of avoiding the use of the measured photon 
energy, which has an asymmetric measurement uncertainty. We improve the measurement 
of the di-hadron 4-momenta (the di-hadron invariant mass resolution becomes 3.2, 2.6, and 
2.0MeV/c 2 for the pion, kaon, and proton modes, respectively) by using the constraint in 
Equation [H and then demanding that: 

0.950 < (E h+h - + £ 7 )/2£ bcam < 1.025. 

Strong electron and muon vetoes are imposed to suppress the abundant QED processes 
e + e~ — > r ye + e~ and e + e~ — > r yfi + fi~. To reject e + e~ — > 7e + e~, we require each track to 
have a matched CC shower with an energy E, together with a measured momentum p, 
such that \E/p — 0.95| > 0.1, and that the combined RICH and dE/dx likelihood for h 
be higher than the combined likelihood for e. To reject e + e~ — > 7// + //~, we require that 
neither track produce a signal in the five hadronic interaction lengths of the muon system. 
For the 7r + 7r~ mode, where muon background is a particular problem because of the similar 
pion and muon masses, we further require that both tracks must be within the barrel part 
of the muon chambers (| cos#| < 0.7), and both have p > 1 GeV/c. To increase the solid- 
angle acceptance of the detector and improve the overall muon suppression efficiency with 
virtually no increase in muon fakes, we flag an event as "not muonic" and remove the muon 
suppression requirements if either track deposits more than 600 MeV in the CC. 

Events that satisfy all the above requirements are then identified as either 7r + 7r _ , K + K~, 
or pp using the RICH and dE/dx information. Since the ratios 7r + 7r~ /K + K~ and K + K~ / pp 
are much larger than 1 for these types of events, in the 3 cases in which we try to reduce 
the background from a lower-mass hadron, we also use the chi-squared value from the kine- 
matic constraint in Equation ^ to identify the event type. Since the constraint involves the 
di-hadron energy, the chi-squared value is sensitive to the hadronic masses. After these pro- 
cedures, the particle identification efficiencies (fake rates) are 90% (0.31%), 99% (0.03%), 
98% (0.10%) for kaons (pions faking kaons), protons (pions faking protons), and protons 
(kaons faking protons), respectively. 

IV. DETERMINATION OF SIGNALS AND THEIR SPIN ASSIGNMENTS 

The overall reconstruction efficiencies as determined by Monte Carlo simulations, includ- 
ing both event selection and analysis cuts, are 43%, 48%, and 56% for the T(1S) radiative 
decays to 7r + 7r~, K + K~ , and pp, respectively. These efficiencies are only mildly dependent 
on the di-hadron invariant mass and are very similar for the continuum background events. 
The continuum-subtracted di-hadron invariant mass plots are obtained by efficiency cor- 
recting each bin of the di-hadron invariant mass plots for the T(1S) and y/s = 10.56 GeV 
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datasets, scaling the latter plot by a factor of 0.404 ± 0.002 4 , and subtracting it from the 
T(1S) dataset invariant mass plot. Possible signals are determined by fitting each spectrum 
to spin- dependent relativistic Breit-Wigner functions 5 . The spin value for each Breit-Wigner 
is surmised by identifying each possible resonance in the invariant mass plot based on its 
approximate mass and width. Later, we confirm these spin assignments for the significant 
resonances by inspecting the angular distributions of the T(1S) decay products. 

The 7r + 7r~ invariant mass plots for the T(1S) and the scaled y/s = 10.56 GeV datasets 
are shown in Figure H The fit to the continuum-subtracted 7r + 7r _ spectrum, shown in 
Figure El has a significant / 2 (1270) signal of 944 ± 74 events. It also has two less significant 
signal candidates; 340l}gQ events in the /o(980) region, and 80 ± 30 events in the ^4(2050) 
region (see Figure EJ) whose significances are 4.3a and 2.6a, respectively. Each significance 
is obtained by doing multiple chi-squared fits to the invariant mass plot fixing the signal 
area to different values, assigning each of these multiple fits a probability proportional to 
e~ x / 2 , normalizing the resulting probability distribution, and calculating the probability for 
negative or signal. 

The K + K~ invariant mass plots for the T(1S) and the scaled y/s = 10.56 GeV datasets 
are shown in Figure 0] The fit to the continuum-subtracted K + K~ spectrum, shown in 
Figure has a significant signal of 312tgi events identified as the / 2 (1525), and two non- 
significant signal candidates indicating possible / 2 (1270) and / (1710) production with 109± 
36 and 73 ± 29 events whose significances are 3.2a and 3.3a, respectively. The excess of 
events in the /2(1270) region is consistent with that expected using the 77r + 7r - data and the 
known branching ratios for the /2(1270). We also note that there is a significant excess of 
220 ± 20 events above 2.0 GeV/c 2 in the K + K~ invariant mass distribution which is not 
associated with any resonant structure. 

The pp invariant mass plots for the T(1S) and the scaled \/s = 10.56 GeV datasets 
are shown in Figure El No recognizable structure is seen in the continuum-subtracted pp 



4 We obtain this factor, /, from the integrated luminosities of the T(1S) and y/s = 10.56 GeV datasets, 
and the assumption that, to first order, the cross sections of the continuum processes in each run are 
proportional to 1/s. This factor is roughly equal to the factor obtained by using the average energy of 
each dataset, 



/ = 0.404 ; 



1.13 fb- 1 /!0.56GeV x2 



3.49 ft.- 1 V 9.46 GeV 



5 The spin-dependent relativistic Breit-Wigner parameterization used has the following probability distri- 
bution for a particular h + h~ invariant mass x > xq, 



dP( x ) oc xx rn T{x) 



where 

r(x) = r 



\ 2S + 1 r,/ N2 

X - X \ 2(x m - x y 



x m -x J (x - x ) 2 + (x m - x ) 2 ' 
In the above expression, x m and r represent respectively the most likely mass and width, and are allowed 
to float during the fit. The values of xq and S are fixed during the fit to the invariant mass threshold 
for the particular mode and the spin of the resonance, respectively. The number of events for each fitted 
signal candidate is obtained by integrating this Breit-Wigner parameterization between threshold and 
3 GeV/c 2 . 
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FIG. 1: Invariant mass of tt + it from e + e — > 77r + 7r for the scaled -^/s = 10.56 GeV dataset (solid 
line), and the T(1S) dataset (circles). The large number of events near 770MeV/c 2 is due to the 
abundant process e + e~ — > jp. 

spectrum, which is shown in Figure In particular, we do not see an enhancement near 
threshold, as might be expected from the BES X(1860) results |57||. There is a non-significant 
excess of 85 ± 18 events in the 2 — 3 GeV/c 2 invariant mass region. 

To confirm the spins of our ^(1270) — > 7r + 7r~ and / 2 (1525) — ► K + K~ signals, we examine 
the absolute value of the cosine of the polar angle of the photon with respect to the beam axis, 
| cos# 7 |, and the absolute value of the cosine of the angle formed by the 3-momentum vector 
of one of the hadrons measured in the di-hadron rest frame with the photon's direction, 
| cos 8h | ■ The event selection efficiency is slightly dependent on both angles, so to minimize 
systematic effects, the | cos 7 1 and | cos 9h | efficiency-corrected distributions are obtained by 
projecting the 2-dimensional bin- by-bin efficiency-corrected (| cos# 7 |, | cos^l) distribution. 
We also subtract the background contributions from the tails of nearby resonances. The 
resulting angular distributions (shown in Figures H and EJ) are simultaneously fit to the 
helicity formalism prediction [H^, |6?1 [73| for different resonance spin hypotheses up to J = 4. 
For the / 2 (1270) the different fit confidence levels are 8 x 1CT 19 , 2 x 1CT 19 , 0.05, 8 x 10~ 12 , 
and 1 x 10~ 12 for the hypotheses J = 0, 1, 2, 3, 4, respectively. For the / 2 (1525) the different 
fit confidence levels are 2 x 10" 4 , 2 x 10~ 4 , 0.23, 8 x 10~ 3 , and 2 x 10~ 3 for the hypotheses 
J = 0, 1,2,3,4, respectively. These results confirm our identification of the resonances, as 



7 




FIG. 2: Invariant mass of tt + tt from T(1S) — > 77r + 7r and the fit to the three spin-dependent 
relativistic Breit-Wigner functions described in the text. 
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FIG. 3: Invariant mass of ir + ir from T(1S) — > 77r + 7r in the region 1.5-3.0 GeV/c 2 . This fit shows 
the small, non-significant, excess found in the /4(2050) region. 
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FIG. 4: Invariant mass of K + K from e + e — * yK + K for the scaled y/s = 10.56 GeV dataset 
(solid line), and the T(1S) dataset (circles). The large number of events near 1.050 GeV/c 2 is due 
to the abundant process e + e~ — ► -yip. 



in both cases the angular distributions of the data strongly favor the J = 2 hypothesis' 
dP0h ^ -|a | 2 x ^(3cos 2 ^- l) 2 x ^(l + cos 2 fl 7 ) + 
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d cos Ohd cos # 7 



15 3 

lad 2 x — sm 2 26 h x -sin 2 6» + ( 2 ) 

16 4 7 

I a, 2 1 2 x — sin 4 Oh x -(1 + cos 2 9. 



16 



T/' 



where a\, A = 0, 1, 2, are the normalized helicity amplitudes, J dPe h ,e 1 — \ a o\ 2 + | fl i| 2 + 
| a-2 1 2 = 1- In other words, \a\\ 2 is the probability of X in T(1S) — > -yX to have helicity ±A. 
Because of the normalization condition, the (9h, 0-y) probability distribution can be described 
by two free parameters, traditionally chosen to be the helicity production ratios, 

\n I 2 \n I 2 

2 °1 j 2 fl 2 

x = - — — and y = - — —. 

\a r \a r 



6 Some authors use a probability distribution that also depends on a third angle, tph [sof . However, extreme 
care must be taken when using this angle because it makes the probability distribution sensitive to the 
relative phases of the helicity amplitudes. Thus, two new free parameters need to be introduced in such 
a probability distribution, as was noted by 0] and correctly implemented by 0, • Otherwise^ the 
measurement of the helicity amplitudes rests on the assumption that their relative phases are |46|, l48l 

El El 
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FIG. 5: Invariant mass of K + K~ from T(1S) — ► r yK + K~ and the fit to the three spin-dependent 
relativistic Breit-Wigner functions described in the text. The dotted line shows the extrapolation 
of this fit to masses above 2 GeV/c 2 , which is the cut-off for the fit. 
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FIG. 6: Invariant mass of pp from e + e — * jpp for the scaled y/s = 10.56 GeV dataset (solid line), 
and the T(1S) dataset (circles). The events near 3.1 GeV/c 2 are due to the process e + e~ — ► r fJ/tp. 
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FIG. 7: Invariant mass of pp from T(1S) — > jpp. 



To measure x 2 and y 2 , we simultaneously fit the data to the individual 9h and # 7 distributions 
using 7 



dN 



d cos 6~ 



dN ff 



d cos 6h 



N / dP t 

N 

1 + x 2 + y 2 

N I dP 9hA 
Je-y 

N 

1 + x 2 + y 2 



;i+y 2 )(l + cos 2 ^ 7 ) + ^ 2 sin 2 ^ 



i IS 

-(3 cos 2 9 h - l) 2 + —x 2 sin 2 + — y 2 sin 4 ^ 
16 " 



15 
— 1 
16' 



(3) 



where N corresponds to the number of events. Using the fits to the data (see Figures |H1 
and EJ) we measure the following helicity production ratios: 



,- _ n nn +0 - 02+0 - 01 n, 2 - n nq+°- 08 + - 04 
x f „fTr,7n\ — u.uu_ 00 _o oo i// 2 (i270) — u.uy_ 



„2 

'/a(1270) 



-0.07-0.03 



™2 n nf) +0.10+0.01 2 _ n on+0.22+0.07 

X /2(1525) _ U - UU -0.00-0.00 i / /2(1525) _ U - OU -0.17-0.06 

where the first uncertainty is statistical and the second is systematic. The systematic un- 
certainty is quantified by studying how well input values are reproduced when analyzing 
Monte Carlo samples, the measured angular distribution of the photon and tracks from 



7 We choose to use a simultaneous fit to these two distributions instead of a two-dimensional fit using 
Equation |21 because of our limited statistics. 
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FIG. 8: Distributions of |cos#,r| (top) and |cos# 7 | (bottom) for the signal events in the /2(1270) 
invariant mass region. The solid lines correspond to a simultaneous fit to the J = 2 helicity 
formalism prediction (Equation |3J) . 



e + e~ — > 7P, 70 in the \/s = 10.56 dataset, the effect of a one-sigma variation in the contin- 
uum scale factor, and possible interference with nearby resonances. We find that possible 
interference with nearby resonances in the | cos 0^1 distribution dominates the systematic 
uncertainty. The helicity production ratio measurements indicate that both resonances are 
predominantly produced with helicity 0. They are in agreement with the predictions of [2], 
and in good agreement with the twist-two-order predictions of yj]: no A = 1 production, and 
A = 2 production suppressed by a factor of (rnx/fn b ) 2 with respect to A = production, 
where mx is the mass of the tensor meson and mj is the mass of the b quark. 

We use the results from fitting the angular distributions to correct the Monte Carlo 
simulation efficiencies, which are calculated using flat distributions in the relevant angles, 
by a factor of 0.78 ± 0.02 for the / a (1270), 0.90 ± 0.01 for the /a (1525), and 0.88^ for 
the significant excess in the 2 — 3GeV/c 2 region of the di-kaon invariant mass. The large 
correction in the pion mode is due to the necessarily stronger muon suppression requirement. 
The measured branching ratios of the significant resonances are: 

B(T(1S) -> 7/2(1270)) = (10.2 ± 0.8 ± 0.7) x 10~ 5 , 
B(T(1S) -> 7/2(1525)) = (3.7l°; 9 7 ± 0-8) x 10~ 5 , 

and the measured branching ratio of the excess events in T(1S) — > , yK + K~ with di-kaon 
invariant mass between 2 — 3 GeV/c 2 is 

B(T(1S) -> -fK + K-) = (1.14 ± 0.08 ± 0.10) x 10~ 5 , 

where the first uncertainty is statistical and the second is systematic. The sources of sys- 
tematic uncertainty are 1% from the number of T(1S) decays, 2% from the Monte Carlo 
simulation of the track reconstruction, 3% (8%) from the Monte Carlo efficiency modeling 
of the event requirements in the pion (kaon) mode, and 1% to 3% from the uncertainty in 
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FIG. 9: Distributions of | cos0r-| (top) and | cos# 7 | (bottom) for the signal events in the / 2 (1525) 
invariant mass region. The solid lines correspond to a simultaneous fit to the J = 2 helicity 
formalism prediction (Equation |3J) . 



the angular distribution measurements. We also assign a 15% systematic uncertainty to the 
/ 2 (1525) branching fraction from possible interference between the /2(1270) and /2(1525) 
resonances, and less than a 1% systematic uncertainty from high- momentum neutral pions 
faking photons in the decay T(1S) — > pn, based on the upper limit in |2f|. Finally, we 
include the uncertainties in the / 2 (1270) and /2(1525) hadronic branching ratios |H| in the 
systematic uncertainty. For our less significant signal candidates, the branching fraction cen- 
tral values, along with their significances and 90% confidence level upper limits, are shown 
in Table |U 



TABLE I: Branching fraction central value (BF), its statistical significance, and its 90% confidence 
level upper limit (UL), for each signal candidate with a significance < 5a. In the branching fraction 
central value, the first uncertainty is statistical and the second is systematic. The first three table 
entries are product branching fractions. 



Channel 


BF (10~ 5 ) 


Significance 


UL (10- 5 ) 


T(1S) - 7 /o(980) -» iT+iT- 
T(1S) -» 7/4(2050) -> ir+ir- 


1.8^ ±0.1 
0.37 ±0.14 ±0.03 


4.3cj 
2.6a 


< 3 
< 0.6 


T(1S) -> 7/o(1710) -» K+K- 


0.38 ±0.16 ±0.04 


3.2a 


< 0.7 


T(1SH IPP, 2GeV/c 2 < m pP <3GeV/c 2 


0.41 ±0.08 ±0.10 


4.8cj 


< 0.6 
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V. DETERMINATION OF UPPER LIMITS FOR /j(2220) AND Z(1860) PRODUC- 
TION AND DECAY 

To measure upper limits of the product branching ratio for the decays T(1S)— > 7/j(2220) 
with /j(2220) — > h + h~, we fit the h + h~ invariant mass plots, shown in Figure ITUl using 
a Breit-Wigner with a peak mass and width fixed at 2.234 GeV/c 2 and 0.017 GeV/c 2 , re- 
spectively. These are the values from the possible /j(2220) signal reported by the BES 
experiment [54|, which is considered a candidate for a glueball. To model the general ex- 
cess of events between 2.0 and 2.5 GeV/c 2 we also use a flat background function in the fit. 
Although the highest bin in the pp plot is indeed in the region of the /j(2220), the excess 
(12 ±5 events) is not significant, and there are no significant signals anywhere in these three 
plots. To find upper limits for /j(2220) — > h + h~ decays, we fix the area of the Breit-Wigner 
to different values, minimize the chi-squared from the fit, and give that area a probability 
proportional e~ x I 2 . These probability distributions are then used to obtain the following 
90% confidence level upper limits on the product branching ratio for /j(2220) production 
and decay to each mode: 

£(T(1S) -> 7 /j(2200)) x £(/j(2200) -> tt+tt-) < 8 x 10~ 7 , 
B(T(1S) -> 7 /j(2200)) x £(/j(2200) -> K + K~) < 6 x 10~ 7 , 
£?(T(1S) - 7 /j(2200)) x i3(/j(2200) - pp) < 11 x lO" 7 . 

The systematic uncertainties on the branching ratios were added in quadrature with the 
statistical errors in forming the above limits. Using the X(1860) parameters measured 
in j5?J, and proceeding in a similar manner as described above, we obtain: 

<B(T(1S) -> 7^(1860) ) x £(X(1860) -> pp) < 5 x 10~ 7 , 

at the 90% confidence level. 

VI. SUMMARY AND CONCLUSION 

We have confirmed CLEO's previous observation of the /2(1270) in radiative T(1S) decays 
and made a new observation of the / 2 (1525), obtaining factors of 0.07±0.01 and 0.081q!o3 f° r 
the ratio of the T(1S) branching fraction with respect to the one measured in J /if} radiative 
decays, respectively. These values are larger than, but the same order of magnitude as, the 
ratio of 0.04 expected from naive scaling arguments. The observed /2(1270) production is in 
agreement with the prediction in and somewhat lower than the prediction in jH . In both 
of the measured modes we can confirm by fits to the angular distributions of the photon and 
charged particles that the two daughter hadrons are indeed produced by a spin-2 parent. 
We find that this parent is produced mostly with helicity 0, in good agreement with the 
predictions in 0, HI . No structure is seen in the pp invariant-mass distribution. In particular, 
we do not observe a near-threshold enhancement as in [5?}. Finally, stringent limits have 
been put on the production of the glueball candidate /j(2220) in radiative T(1S) decays. 
Glueball production is expected to be enhanced in T(1S) radiative decays 0,173, [lH, but 
we find that, within our experimental sensitivity, known tensor meson states, believed to be 
composed only of quarks, dominate the di-gluon spectrum. 
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FIG. 10: Invariant mass of ir + ir~ (top) K + K~ (middle) and pp (bottom) from T(1S) — »■ jh + h~ in 
the range 2.0 to 2.5 GeV/c 2 . Also shown are the results of fits to a /j(2220) Breit-Wigner function 
(dashed) and a flat background (solid). The /j(2220) signal line shown corresponds to the 90% 
confidence level maximum yield obtained from the fit. The flat background function is also used 
to measure the possible excesses of events in the invariant mass region 2 to 3 GeV/c 2 . 
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